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PREFACE 


The  investigation  reported  herein  was  authorized  by  the  Office, 
Chief  of  Engineers  (OCE),  U.  S.  Army,  in  a letter  dated  2U  June  1975 
and  was  performed  for  the  Federal  Insurance  Administration,  Department 
of  Housing  and  Urban  Development,  under  Inter-Agency  Agreements  IAA-H- 
16-75,  Project  Order  No.  23;  IAA-H-7-76,  Project  Order  No.  2;  and  IAA 
H-10-77 , Project  Order  No.  30.  Project  coordinator  was  Mr.  Jerome 
Peterson,  OCE. 

The  investigation  was  conducted  from  July  1975  to  November  1977 
by  personnel  of  the  Hydraulics  Laboratory,  U.  S.  Army  Engineer  Water- 
ways Experiment  Station  (WES),  under  the  direction  of  Mr.  H.  B.  Simmons, 
Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  Chief  of  the  Wave 
Dynamics  Division,  and  Mr.  D.  D.  Davidson,  Chief  of  the  Wave  Research 
Branch.  Messrs.  J.  R.  Houston,  Research  Physicist,  and  A.  W.  Garcia, 
Research  Oceanographer,  both  of  the  Hydraulics  Laboratory,  conducted 
the  study;  and  this  report  was  prepared  by  Mr.  Houston. 

Directors  of  WES  during  the  investigation  and  the  preparation 
and  publication  of  this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  John  L. 
Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI)  AND 
METRIC  (SI)  TO  U.  S.  CUSTOMARY  UNITS  OF  MEASUREMENTS 


P 


Units  of  measurement  used  in  this  report  can  be  converted  as  follows: 


Multiply 


To  Obtain 


U.  S.  Customary  to  Metric  (Si) 


feet  0.30U8 

feet  per  second  per  second  0.30U8 

miles  (U.  S.  statute)  1.852 

miles  (U.  S.  nautical)  1.852 

per  hour 

square  feet  per  second  0.09290301* 

per  second 


metres 

metres  per  second  per  second 

kilometres 

kilometres  per  hour 

square  metres  per  second 
per  second 


kilometres 

metres 


Metres  (Si)  to  U.  S.  Customary 

0.6213711  miles  (U.  S.  statute) 

3.280839  feet 
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TYPE  16  FLOOD  INSURANCE  STUDY:  TSUNAMI  PREDICTIONS 


FOR  THE  WEST  COAST  OF  THE  CONTINENTAL 
UNITED  STATES 


PART  I : INTRODUCTION 

Background 


1.  Of  all  water  waves  that  occur  in  nature,  one  of  the  most  de- 
structive is  the  tsunami.  The  term  "tsunami,"  originating  from  the 
Japanese  words  "tsu"  (harbor)  and  "nami"  (wave),  is  used  to  describe 
sea  waves  of  seismic  origin.  Tectonic  earthquakes,  i.e.,  earthquakes 
that  cause  a deformation  of  the  seabed,  appear  to  be  the  principal  seis- 
mic mechanism  responsible  for  the  generation  of  tsunamis.  Coastal  and 
submarine  landslides  and  volcanic  eruptions  also  have  triggered 
tsunamis. 

2.  Tsunamis  are  principally  generated  by  undersea  earthquakes  of 
magnitudes  greater  than  6.5  on  the  Richter  scale  with  focal  depths  less 
than  50  km.#  They  are  very  long-period  waves  (5  min  to  several  hours ) 
of  low  height  (a  few  feet  or  less)  when  traversing  water  of  oceanic 
depth.  Consequently,  they  are  not  discernible  in  the  deep  ocean  and 

go  unnoticed  by  ships.  Tsunamis  travel  at  the  shallow-water  wave 
celerity  equal  to  the  square  root  of  acceleration  due  to  gravity  times 
water  depth  even  in  the  deepest  oceans  because  of  their  very  long  wave- 
lengths. This  speed  of  propagation  can  be  in  excess  of  500  mph  in  the 
deep  ocean. 

3.  When  tsunami  waves  approach  a coastal  region  where  the  water 
depth  decreases  rapidly,  wave  refraction,  shoaling,  and  bay  or  harbor 
resonance  may  result  in  significantly  increased  wave  heights.  The  great 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Si)  units  find  metric  (Si)  units  to  U.  S.  customary 
is  presented  on  page  4. 
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period  and  wavelength  of  tsunami  waves  preclude  their  dissipating  energy 
as  a breaking  surf;  instead,  they  are  apt  to  appear  as  rapidly  rising 
water  levels  and  only  occasionally  as  bores. 

U.  Over  500  tsunamis  have  been  reported  within  recorded  history. 
Virtually  all  of  these  tsunamis  have  occurred  in  the  Pacific  Basin. 

This  is  because  most  tsunamis  are  associated  with  earthquakes,  and 
most  seismic  activity  beneath  the  oceans  is  concentrated  in  the  narrow 
fault  zones  adjacent  to  the  great  oceanic  trench  systems  which  are  pre- 
dominantly confined  to  the  Pacific  Ocean. 

5.  The  loss  of  life  and  destruction  of  property  due  to  tsunamis 
have  been  immense.  The  Great  Hoei  Tokaido-Nankaido  tsunami  of  Japan 
killed  30,000  people  in  1707-  In  1868,  the  Great  Peru  tsunami  caused 
25,000  deaths  and  carried  the  frigate  U.  S.  S.  Waterlee  1,300  ft  inland. 
The  Great  Meiji  Sanriku  tsunami  of  1896  killed  27,122  persons  in  Japan 
and  washed  away  over  10,000  houses. 

6.  In  recent  times,  three  tsunamis  have  caused  major  destruction 
in  areas  of  the  United  States.  The  Great  Aleutian  tsunami  of  19^6 
killed  173  persons  in  Hawaii,  where  heights  as  great  as  55  ft  were 
recorded.  The  i960  Chilean  tsunami  killed  330  people  in  Chile,  6l  in 
Hawaii,  and  199  in  distant  Japan.  The  most  recent  major  tsunami  to 
affect  the  United  States,  the  196k  Alaskan  tsunami,  killed  107  people 
in  Alaska,  U in  Oregan,  and  11  in  Crescent  City,  California,  and  caused 
over  100  million  dollars  in  damage  on  the  west  coast  of  North  America. 

Purpose  of  Study 

7.  The  purpose  of  this  study  was  to  establish  100-  and  500-year 
tsunami  runup  elevations  on  the  west  coast  of  the  continental  United 

1 2 

States  produced  by  distantly  generated  tsunamis . Two  earlier  reports  ’ 
established  these  runup  elevations  for  Southern  California,  Monterey 
Bay,  San  Francisco  Bay,  and  Puget  Sound.  The  100-  and  500-year  tsunami 
runup  elevations  are  required  by  the  Federal  Insurance  Administration 
(FIA)  of  the  Department  of  Housing  and  Urban  Development  for  use  in 
flood  insurance  rate  calculations. 
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PART  II:  EXPLANATION  OF  RESULTS 


8.  Plates  1-30  present  predicted  100-  and  500-year  elevations 
(in  feet)  produced  by  distantly  generated  tsunamis  on  the  west  coast  of 


the  continental  United  States.  These  elevations  include  the  effects 


of  the  astronomical  tide;  that  i3,  they  are  maximum  elevations  due  to 
the  superposition  of  tsunami  and  tidal  wave  forms  (see  PART  IV).  The 


lower  curves  in  Plates  1-30  represent  the  100-year  runup  and  the  upper 


curves,  the  500-year  runup.  A 100-year  runup  is  one  that  is  equaled  or 
exceeded  with  an  average  frequency  of  once  every  100  years ; a 500-year 
runup  has  a corresponding  definition.  Runup  values  in  this  report  are 
referenced  to  the  mean  sea  level  (msl)  datum. 

9.  Plates  1 and  2 present  100-  and  500-year  runup  elevations 
versus  longitude  in  minutes.  These  two  plates  only  cover  a section  of 
coast  extending  from  longitude  119°45'  to  120°30',  (from  Santa  Barbara, 
California,  to  latitude  34°30').  This  is  a section  of  coast  oriented 


east-west.  The  remainder  of  the  west  coast  has  a north-south  orienta- 


tion and  100-  and  500-year  runup  elevations  versus  latitudes  in  minutes 
for  this  coast  (extending  from  3/*°30'  to  48°23.6')  are  presented  in 


Plates  3-30. 


10.  In  order  to  determine  100-  or  500-year  runup  elevations  at  a 


coastal  location,  it  is  necessary  to  know  the  latitude  of  the  location 
(or  the  longitude  of  the  location  if  the  latitude  is  less  than  31*°30'). 
The  latitude  and  longitude  of  a coastal  location  can  be  determined 
using  most  standard  maps.  For  example,  the  National  Ocean  Survey 
(previously  the  U.  S.  Coast  and  Geodetic  Survey)  maps  and  topographic 
quadrangle  maps  published  by  the  U.  S.  Geological  Survey  have  latitude 
and  longitude  markings. 

11.  Figure  1 shows  a section  of  coast  in  northern  California 
taken  from  the  U.  S.  Geological  Survey  quadrangle  map  entitled  "Cannibal 
Island,  California."  This  quadrangle  extends  from  a latitude  of 
1*0°3T'30"  to  U0°U5’.  The  section  of  coast  shown  in  Figure  1 extends 
from  U0°U0'  to  U0°U2'30"  with  the  40'  and  42' 30"  markings  located  on 
the  right-hand  side  of  the  figure.  The  40°  marking  is  on  a section  of 
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the  quadrangle  not  shown.  If  the  runup  elevations  are  needed  at  the 
Table  Bluff  Lighthouse  shown  in  Figure  1,  the  latitude  of  this  coast 
location  must  first  be  determined.  Using  an  ordinary  ruler  and  the 
latitude  markings  given,  the  latitude  of  the  Table  Bluff  Lighthouse  is 
found  to  be  l*0°Ul.7".  The  runup  elevations  for  this  latitude  are  pre- 
sented in  Plate  15.  The  100-year  runup  elevation  is  10.2  ft  and  the 
500-year  runup  elevation  is  20.6  ft. 

12.  Several  bays  along  the  west  coast  have  multiple-valued  coast- 
lines (that  is,  within  and  near  a bay  there  are  locations  having  the 
same  latitude);  these  are  San  Luis  Obispo,  Bolinas,  Drake's,  Bodega, 
Trinidad,  and  Crescent  City  Bays  in  California  and  Port  Orford  Bay  in 
Oregon.  Figure  2 is  an  illustrative  drawing  of  such  bays.  Note,  for 
example,  that  locations  h,  6,  and  9 have  the  same  latitude.  Thus,  the 
plots  of  runup  elevation  versus  latitude  must  necessarily  be  multiple- 
valued. Figure  2 shows  that  locations  h,  6,  and  9 have  the  same  lati- 
tude but  different  elevation  values.  To  eliminate  confusion  as  to  which 
runup  elevation  corresponds  to  a particular  latitude,  the  following 
convention  is  used:  movement  along  a coast  from  south  to  north 
corresponds  to  movement  from  left  to  right  along  the  runup  curves  in 
Plates  3-30.  Figure  2 illustrates  this  convention  using  numbers  to 
represent  locations. 

13.  This  report  does  not  present  runup  elevations  for  the  Cali- 
fornia coast  between  San  Diego  and  Santa  Barbara  since  such  elevations 
were  presented  in  an  earlier  report.^-  Runup  elevations  for  Monterey 
Bay,  San  Francisco  Bay,  and  Puget  Sound  were  also  presented  in  an 
earlier  report.2 


LATITUDE 


Figure  2.  Multi pie- valued  coastlines 


Background 


lU.  Unlike  other  areas  of  the  Pacific  Ocean  such  as  the  Hawaiian 
Islands,  the  west  coast  of  the  continental  United  States  lacks  suffi- 
cient data  to  allow  tsunami  elevation  predictions  based  upon  local 
historical  records  of  tsunami  activity.  Virtually  all  of  the  west 
coast  is  completely  without  data  of  tsunami  occurrence,  even  data  for 
the  prominent  196U  tsunami.  Only  a handful  of  locations  have  histori- 
cal data  for  tsunamis  other  than  the  196U  tsunami.  However,  the  FIA 
requires  information  on  tsunami  elevations  for  the  entire  west  coast 
of  the  continental  United  States,  even  for  the  many  locations  that 
have  no  known  historical  data  of  tsunami  activity  and  for  coastal 
areas  that  are  currently  not  developed  (since  these  areas  may  be 
developed  in  the  future ) . 

15*  The  lack  of  historical  data  of  tsunami  activity  on  the  west 
coast  necessitates  the  use  of  numerical  models  to  predict  runup  eleva- 
tions. The  Aleut ian-Alaskan  area  and  the  west  coast  of  South  America 
were  found  in  Reference  1 to  be  the  tsunamigenic  regions  of  concern  to 
the  west  coast  of  the  United  States.  Both  regions  have  sufficient  data 
on  the  generation  of  major  tsunamis  to  allow  a statistical  investigation 
of  tsunami  generation.  A numerical  model  employing  a fairly  crude 
grid  to  cover  a large  section  of  the  Pacific  Ocean  was  used  to  generate 
representative  tsunamis  and  propagate  them  across  the  deep  ocean.  A 
second  numerical  nv  tlel  employing  a fine  grid  model  was  used  to  propa- 
gate tsunamis  from  the  deep  ocean  over  the  continental  slope  and  shelf 

1 2 

to  shore.  Previous  studies  ’ used  analytical  methods  to  solve  one- 
dimensional linear  equations  that  described  the  transformation  of 
tsunamis  as  they  propagated  from  deep  water  to  shore.  The  two- 
dimensional  numerical  model  used  in  this  study  to  propagate  tsunamis  to 
shore  solved  equations  that  included  nonlinear  advective  terms  and 
bottom  stress  terms. 

16.  In  this  report,  only  tsunamis  of  distant  origin  (at  least  two 


11 


or  three  tsunami  wavelengths  away  from  the  west  coast)  are  considered 

3 

in  the  analysis.  Hammack  has  shown  that  near  the  generation  area  of 

a tsunami,  details  of  ground  motion  during  the  earthquake  and  details 

of  the  permanent  deformation  of  the  seafloor  influence  the  form  of  the 

resulting  tsunami.  Very  little  is  known  about  the  actual  time-dependent 

ground  motion  during  earthquakes  generating  major  tsunamis,  and  small- 

scale  details  of  the  permanent  deformation  of  the  seafloor  following 

earthquakes  cannot  be  predicted  in  advance.  Thus,  accurate  predictions 

of  the  properties  of  locally  generated  tsunamis  are  not  possible  at 

3 

this  time.  However,  Hammack  has  shown  that  the  time-dependent  ground 
motions  and  small-scale  details  of  the  permanent  ground  deformation 
produce  waves  which  are  not  significant  far  from  the  source  region. 

Thus,  distantly  generated  tsunamis  can  be  studied  only  knowing  major 
features  of  the  permanent  ground  deformation. 

17-  The  probability  is  not  considered  very  great  that  a destruc- 
tive, locally  generated  tsunami  will  occur  on  the  west  coast  of  the 
continental  United  States.  Tsunamis  are  generally  produced  by  earth- 
quakes having  fault  movements  that  exhibit  a pronounced  "dip-slip," 
or  vertical  component  of  motion.  "Strike-slip,"  or  horizontal  displace- 
ment, fault  movements  are  inefficient  generators  of  tsunamis.  Faults 
on  the  west  coast  of  the  United  States  characteristically  exhibit 
strike-slip  motion  since  the  Pacific  block  of  the  earth’s  crust  is 
moving  horizontally  relative  to  the  North  American  block.  The  west 
coast  of  the  United  States  does  not  share  the  characteristics  (oceanic 
trenches  and  island  arcs)  of  known  tsunami-generating  areas  and,  in 
fact,  has  not  historically  been  one.  Relatively  small  locally  generated 
tsunamis  have  been  known  to  occur  on  the  west  coast , but  there  are  no 
reliable  reports  of  major  locally  generated  tsunamis.  There  could  be 
a few  locations  on  the  west  coast  for  which  locally  generated  tsunamis 
pose  a greater  hazard  than  do  distantly  generated  tsunamis  because  the 
elevations  produced  by  distantly  generated  tsunamis  are  small.  However, 
predictions  of  elevations  produced  by  locally  generated  tsunamis  are 
beyond  the  scope  of  this  report. 


Generation  and  Deep-Ocean  Numerical  Model 

18.  The  finite-difference  numerical  model  used  to  simulate  the 
generation  of  tsunamis  and  their  propagation  across  the  deep  ocean  was 
originally  developed  by  Hwang  et  al. , and  described  in  detail  in  an 
earlier  report.1  In  this  study  the  model  employed  l/3°  by  1/3°  spheri- 
cal coordinate  grids  to  solve  the  linearized  long-wave  equations.  Such 
crude  grids  adequately  resolve  the  very  long  tsunami  wavelengths  in 
the  deep  ocean.  A comparison  of  a wave  gage  recording  of  the  196^ 
tsunami  at  Wake  Island  in  the  Pacific  Ocean  and  a simulation  of  the 
1964  tsunami  using  this  numerical  model  and  1/3°  by  1/3°  grid  verified 
the  accuracy  of  such  a grid  spacing."^  The  grids  covered  very  large 
sections  of  the  Pacific  Ocean  including  either  Alaska  and  the  west 
coast  of  the  United  States  or  South  America  and  the  west  coast  of  the 
United  States.  The  boundary  condition  on  the  solid  boundaries  (land) 
of  the  grids  was  that  the  component  of  the  velocity  normal  to  the  bound- 
ary equals  zero.  On  open  boundaries  (ocean),  a first-order  approxima- 
tion of  total  transmission  was  made. 

19*  The  deep-ocean  finite-difference  model  solved  an  initial 
value  problem  starting  with  an  uplift  deformation  of  the  water  surface 
identical  with  the  major  features  of  the  permanent  deformation  (perma- 
nent in  the  sense  that  the  time  scale  associated  with  it  is  much  longer 
than  the  period  of  the  tsunami ) of  the  seafloor  following  the  seismic 

disturbance.  The  transient  movements  within  the  time  history  of  the 

3 

ground  motion  werfe  neglected  because  Hammack  has  shown  that  these 
movements  are  unimportant  in  the  far  field  for  a spatially  large  impul- 
sive ground  motion.  Hammack  has  further  shown  that  the  inital  deforma- 
tion of  the  water  surface  will  closely  approximate  major  features  of 
the  permanent  deformation  of  the  ocean  floor,  provided  these  features 
have  characteristic  lengths  that  are  at  least  four  times  as  great  as 
the  water  depth.  The  neglect  of  smaller  features  is  unimportant 
because  such  small-scale  details  produce  waves  that  are  negligible  in 
the  far  field. ^ 
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Nearshore  Numerical  Model 


20.  The  nearshore  finite-difference  numerical  model  was  based 

7 

upon  the  original  formulation  of  a tidal  hydraulics  model  by  Leendertse. 

The  model  has  been  applied  in  several  studies  to  analyze  the  tidal 

8 9 

hydraulics  of  harbors  and  inlets.  * It  also  has  been  used  to  calculate 

2 

decay  of  a tsunami  as  it  enters  and  spreads  throughout  San  Francisco  Bay. 

21.  The  fundamental  equations  solved  by  this  numerics!  model  were 
vertically  arranged  and  expressed  in  a Cartesian  coordinate  system.  The 
momentum  equations  were 


w + „ au  au  ia  + _jb_ 

at  ax  ay  K ax  h + n 


= o 


B 


3V  + U 3V  + V3V  la* 

at  3x  3y  3y  h + n 


= 0 


where 


2 2 
U + V 


T*  = gU/C^ 

= gv/c2>ju2  + v: 

T*  and  T^  = bottom  stress  terms 

O D 

C = Chezy  coefficient 
The  continuity  equation  was: 


(1) 

(2) 


!?  * k <h  * ")D  * k <h  4 n)v  ‘ 0 


(3) 


23.  This  set  of  equations  was  solved  by  a system  of  finite- 
difference  equations  using  central  differences  on  a space-staggered 
grid.  Leendertse' s implicit-explicit  multi operational  method  was 
employed  in  determining  n as  a function  of  time. 

2b.  The  nearshore  numerics!  model  used  as  input  the  time-history 
calculated  by  the  generation  and  deep-ocean  propagation  numerical  model. 
A tsunsuni  was  generated  in  the  Aleutietn-Alaskan  area  or  the  west  coast 
of  South  America  and  propagated  across  the  deep  ocean  to  a 500-m  depth 
off  the  west  coast  of  the  United  States.  Wave  forms  calculated  at  this 


depth  by  the  deep-ocean  numerical  model  were  recorded  all  along  the 
west  coast.  These  wave  forms  then  were  used  as  input  to  the  nearshore 
numerical  model,  which  propagated  the  tsunamis  from  the  500-m  depths 
across  the  continental  slope  and  shelf  to  shore. 

25.  Figure  3 shows  outlines  of  the  four  numerical  grids  used  to 
cover  the  west  coast.  The  grids  had  square  grid  cells  2 miles  on  a 
side.  The  offshore  bathymetry  was  modeled  from  the  500-m  contour  to 
shore.  Beyond  the  500-m  contour,  the  ocean  was  modeled  as  being  a con- 
stant 500-m  depth.  The  input  boundary  of  each  grid  was  located 
approximately  one  and  one-half  wavelengths  of  a 30-min  wave  from  the 
3hore.  Therefore,  at  least  three  typical  waves  could  arrive  at  the 
shore  before  waves  were  rereflected  from  the  input  boundary.  The  wave 
forms  used  as  forcing  functions  varied  along  the  input  boundaries  in 


Figure  3.  Nearshore  grids 
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accordance  with  the  calculations  of  the  generation  and  deep-ocean  prop- 
agation model.  The  input  boundaries  of  the  grids  were  oriented  approx- 
imately parallel  to  the  shoreline,  since  refraction  would  have  bent  the 
wave  fronts  to  such  an  orientation  before  they  reached  the  500-m  con- 
tour.10 Lateral  boundaries  of  the  grids  were  impermeable  vertical 
walls. 

1 2 

26.  In  previous  studies,  ’ an  analytical  solution  was  used  to 
propagate  tsunamis  across  the  continental  slope  and  shelf  to  shore. 

This  analytical  approach  had  limitations  since  it  solved  a one- 
dimensional, linear,  dissipation- free  equation  and  required  the  wave 
form  to  be  sinusoidal  and  of  constant  height.  The  finite-difference 
model  used  in  this  study  remedied  these  limitations  since  it  solved 
two-dimensional  equations  that  included  both  nonlinear  advective  terms 
and  bottom  dissipation  terms  and  it  also  allowed  input  of  arbitrary 
wave  forms . 

27.  A recent  study  predicting  tsunami  runup  elevations  in  the 

11  12 

Hawaiian  Islands  ’ used  a finite-element  numerical  model  for  near- 
shore tsunami  propagation.  However,  this  model  used  a time-harmonic 
solution  of  linear  and  dissipation- free  equations.  Arbitrary  wave 
forms  could  be  considered  by  the  model,  but  only  through  Fourier  decom- 
position of  the  arbitrary  wave  form,  separate  response  calculations  for 
each  component,  and  then  a Fourier  recomposition  of  the  product  of  the 
components  and  the  responses  calculated  by  the  model.  The  finite 
element  approach  was  used  to  study  Hawaiian  tsunamis  because  these 
islands  are  so  small  and  the  bathymetry  surrounding  them  is  so  rapidly 
varying  that  only  a telescoping  finite  element  grid  could  adequately 
represent  important  physical  features  of  the  problem.  It  would  not 
be  feasible  for  the  nearshore  finite  difference  numerical  model  to  use 
a grid  spacing  as  small  as  that  employed  by  the  finite  element  model 
and  still  cover  all  of  the  Hawaiian  Islands.  The  remarkably  narrow 
continental  shelf  of  the  islands  also  limits  the  time  available  for 
nonlinear  and  dissipative  effects  to  become  important.  Thus,  the  finite 
element  model  was  ideal  for  studying  tsunami  interaction  with  the 
Hawaiian  Islands,  but  of  less  value  for  studying  tsunamis  on  the  west 
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coast  of  the  United  States  where  land-water  boundary  changes  and  depth 
variations  are  relatively  gradual  and  nonlinear  and  dissipative  effects 
possibly  of  greater  importance. 

Verification 


28.  The  numerical  models  used  in  this  report  were  verified  by 
numerical  simulations  of  the  1964  Alaskan  tsunami.  This  was  the  only 
large  tsunami  for  which  reliable  information  exists  concerning  source 
characteristics.  The  initial  condition  used  in  the  generation  and 
deep-ocean  propagation  numerical  model  was  that  the  uplift  of  the 
ocean's  surface  in  the  source  region  was  identical  with  the  permanent 
deformation  of  the  ocean  bottom  following  the  earthquake.  The  perma- 
nent deformation  of  the  ocean's  bottom  as  a function  of  spatial  location 

13 

was  taken  from  Plafker. 

29.  Figure  1+  shows  surface  elevation  contours  approximately 
3-1/2  hr  after  the  1964  Alaskan  earthquake,  as  calculated  by  the  deep- 
ocean  propagation  numerical  model.  This  figure  illustrates  the  concen- 
tration of  energy  on  the  northern  California,  Oregon,  and  Washington 
coasts  due  to  the  directional  radiation  of  energy  from  the  source  re- 
gion. The  upper  section  of  Figure  5 shows  a time-history  of  the  1964 
Alaskan  tsunami  calculated  by  the  deep-ocean  numerical  model  near 
Crescent  City,  California,  in  a water  depth  of  500  m.  This  wave  form 
was  used  as  input  to  the  nearshore  propagation  numerical  model. 

30.  The  lower  section  of  Figure  5 shows  a comparison  between  a 
tide  gage  recording  of  the  1964  tsunami  at  Crescent  City  and  the  near- 
shore numerical  model  calculations.  The  numerical  model  calculations 
agree  quite  well  with  the  tide  gage  recording.  Periods,  phases,  and 
amplitudes  are  in  good  agreement.  The  greatest  disagreement  is  the 
amplitude  of  the  third  wave.  It  should  be  noted  that  there  is  disagree- 
ment on  the  amplitude  of  the  actual  historical  third  and  fourth  wave 
crests.  Figure  5 shows  a reconstruction  of  the  1964  tsunami  at  Crescent 

City  inferred  by  Wilson  and  Torum10  from  the  prototype  gage  record, 

14 

and  later  survey  measurements  by  Magoon.  The  elevation  of  the  fourth 
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Figure  4.  Water-surface  elevations  3-1/2  hr 
after  196U  earthquake  (Hwang,  1972) 

wave  (not  shown)  was  estimated  by  Wilson  and  Torum  to  have  been  a lit- 
tle less  than  lU  ft  above  mean  lower  low  water  (mllw).  The  zero  eleva- 
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tion  shown  in  Figure  5 is  the  mllw  datum.  Wiegel,  however,  estimated 
that  the  elevation  of  the  third  wave  was  approximately  16  ft  above  mllw 
and  that  the  fourth  wave  attained  the  highest  elevation  of  18  or  19  ft 
above  mllw  at  the  tide  gage  location.  The  numerical  model  predicted  a 
fourth  wave  elevation  approximately  the  same  as  that  estimated  by  Wil- 
son and  Torum  (ih  ft  above  mllw).  However,  reflections  from  the  input 
boundary  were  probably  growing  in  importance  in  the  numerical  model 
calculations  during  the  arrival  of  this  fourth  wave. 

31.  There  may  be  several,  reasons  why  the  maximum  wave  elevation 
predicted  by  the  numerical  model  is  not  in  as  good  agreement  with  the 
historical  record  as  are  the  elevations  of  the  initial  waves.  First, 
the  later  waves  in  the  deepwater  wave  form  are  probably  waves  which 
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have  been  reflected  from  land  areas  in  the  source  region.  The  fairly 
large  grid  cell  spacing  used  in  the  deep-ocean  propagation  model  to 
represent  the  land-water  boundary  in  this  region  probably  distorts  the 
reflected  wave  form  somewhat.  Total  reflection  of  this  region  also  was 
not  completely  realistic.  Furthermore,  the  nearshore  grid  may  not  have 
fine  enough  grid  cells  to  completely  model  the  resonance  phenomenon 
leading  to  the  large  third  or  fourth  wave. 

32.  Since  196U,  there  has  been  considerable  speculation  concern- 
ing the  reasons  that  the  effects  of  the  196U  tsunami  were  so  great  at 
Crescent  City.  The  finite-difference  numerical  models  showed  that  both 
the  directional  radiation  of  energy  from  the  source  region  (Figure  4) 
and  a local  resonance  caused  the  relatively  large  elevations  at  Crescent 
City.  Actually,  the  large  elevations  at  Crescent  City  were  not  unique 
since  directional  radiation  of  energy  from  the  source  also  caused  large 
elevations  along  the  Oregon  and  northern  California  coasts.  Runup  10 

to  15  ft  above  the  high  tide  level  occurred  all  along  the  Oregon  coast 
south  of  the  Columbia  River  (runup  at  Crescent  City  was  approximately 
15  ft  above  the  high  tide  le'el).  However,  the  Oregon  coast  is  very 
sparsely  populated  and  thus  there  were  few  damage  reports.  The  severity 
of  structural  damage  at  Crescent  City,  which  has  a large  logging  indus- 
try, also  was  apparently  due  to  the  impact  of  logs  carried  by  the 
tsunami.10  The  nearshore  numerical  model  indicated  that  the  resonant 
effects  at  Crescent  City  were  fairly  local,  extending  over  2 to  U miles 
of  coastline.  This  is  in  disagreement  with  Wilson  and  Torum's10  specu- 
lation that  a bowl-shaped  section  of  the  continental  shelf  with  a diam- 
eter of  approximately  50  miles  experienced  a resonant  oscillation.  His- 
torical data  support  the  numerical  model  calculations.  At  the  mouth  of 
the  Klamath  River,  approximately  15  miles  south  of  Crescent  City,  ele- 
vations observed  during  the  196U  tsunami  were  only  2 to  3 ft  above  nor- 
mal high  tide.  Wilson  and  Torum’s  hypothesis  of  a shelf  oscillation 
would  predict  elevations  at  the  mouth  of  the  Klamath  River  greater  than 
those  that  occurred  at  Crescent  City. 

33.  Figure  6 shows  a comparison  between  a tide  gage  recording  of 
the  196U  tsunami  at  Avila  Beach,  California,  and  the  nearshore  model 
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Figure  6.  196U  tsunami  from  Alaska  recorded 

at  Avila  Beach,  California 

calculations.  The  elevations  recorded  by  the  Avila  Beach  tide  gage 
were  larger  than  those  recorded  at  any  tide  gage  on  the  west  coast 
except  the  Crescent  City  gage.  The  historical  and  calculated  wave  forms 
shown  in  Figure  6 are  in  good  general  agreement.  The  tide  gage  record 
obviously  has  higher  frequency  components  not  predicted  by  the  numeri- 
cal model.  These  components  may  be  local  oscillations  of  water  areas 
which  are  too  small  to  be  accurately  represented  by  the  numerical  grid. 
Important  features  such  as  the  wave  amplitudes  are  in  good  agreement. 

3U.  The  good  agreement  between  the  numerical  model  simulations 
and  tide  gage  recordings  of  the  196U  tsunami  at  the  two  locations  where 
the  largest  waves  were  recorded  in  196^  by  tide  gages  verified  the 
numerical  techniques  described  in  the  previous  section.  Other  tide  gage 
recordings  of  the  196U  tsunami  on  the  west  coast  (except  at  Astoria, 
Oregon)  were  in  areas  not  covered  by  the  grids  shown  in  Figure  3 (these 
areas  were  considered  in  References  1 and  2).  A comparison  with  the 
tide  gage  at  Astoria,  Oregon,  was  not  made  because  Astoria  is  approxi- 
mately 12  miles  away  from  the  coast  in  the  estuary  of  the  Columbia 
River.  Waves  of  the  1961*  tsunami  were  small  at  Astoria. 
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PART  IV:  METHODOLOGY  FOR  RUNUP  PREDICTIONS 


r 


Tsunami  Occurrence  Probabilities 

35.  Historical  data  on  tsunami  generation  must  be  the  basis  for 
an  analysis  that  considers  the  probability  of  tsunami  generation  in 
the  two  tsunamigenic  areas  in  the  Pacific  Ocean  of  concern  to  the  west 
coast  of  the  continental  United  States — the  Aleutian  and  Peru-Chile 
Trench  areas.  A satisfactory  correlation  between  earthquake  magnitude 
and  tsunami  intensity  has  never  been  demonstrated.  Not  all  large  earth- 
quakes occurring  in  the  ocean  even  generate  noticeable  tsunamis.  Fur- 
thermore, earthquake  parameters  of  importance  to  tsunami  generation, 
such  as  focal  depth  and  vertical  ground  motion,  have  only  been  measured 
for  earthquakes  occurring  in  recent  years.  Therefore,  data  on  earth- 
quake occurrence  cannot  be  used  to  determine  occurrence  probabilities 

of  tsunamis.  Historical  data  of  tsunami  occurrence  in  generation  re- 
gions must  be  used  to  determine  these  probabilities. 

36.  In  South  America,  a wealth  of  information  exists  concerning 
tsunami  occurrence.  Reliable  data  (grouped  in  intensity  increments  of 
0.5)  existed  for  tsunamis  with  intensity  greater  than  or  equal  to  0 for 
a 169-vear  period  and  greater  than  or  equal  to  2.5  for  a Ul7-year 
period.  The  intensity  scale  used  was  a modification  (by  S.  L.  Soloviev 
of  the  Soviet  Union)  of  the  standard  Imamura-Iida  tsunami  intensity 
scale.  Intensity  is  defined  as 

i - log2  HavB) 

This  definition  in  terms  of  an  average  runup  (in  metres)  over  a coast 
instead  of  a maximum  runup  elevation  at  a single  location  (used  for  the 
standard  Imamura-Iida  scale)  tends  to  eliminate  a spurious  intensity 
magnitude  caused  by  often  observed  anomalous  responses  (due,  for  exam- 
ple, to  local  resonances)  of  single  isolated  locations. 

37.  Using  the  most  recent  and  complete  catalog  of  tsunami  occur- 
rence in  the  Pacific  Ocean, ^ a relation  between  tsunami  intensity 
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and  frequency  of  occurrence  was  determined  for  the  tsunami -generating 
trench  running  the  length  of  the  Peru-Chile  coast.  Tsunamis  with  in- 
tensity greater  than  or  equal  to  0 were  considered.  It  was  assumed 
that  the  logarithm  of  the  tsunami  frequency  of  occurrence  was  linearly 

related  to  the  tsunami  intensity.  Earthquake  magnitude  and  frequency 
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of  occurrence  have  been  similarly  related  by  Gutenberg  and  Richter 

1 O 

and  used  extensively  in  earthquake  predictions.  Soloviev"^  has  shown 
a similar  relation  between  tsunami  intensity  and  frequency  of  genera- 
tion for  moderate  to  large  tsunamis  throughout  the  Pacific  Ocean.  Fur- 
thermore, Wiegel1^  found  the  same  type  of  relation  between  tsunami 

occurrence  and  runup  levels  for  historical  tsunamis  at  Hilo,  Hawaii; 
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San  Francisco,  California;  and  Crescent  City,  California;  and  Adams 

for  tsunamis  at  Kahuku  Point,  Oahu.  A recent  study  by  Rascon  and 
20 

Villarreal  revealed  this  same  relation  for  historical  tsunamis  on  the 
west  coast  of  Mexico  (data  from  1732)  and  on  the  Pacific  west  coast  of 
North  America,  excluding  Mexico. 

38.  Letting  n(i)  equal  the  probability  of  a tsunami  with  an 
intensity  i being  generated  during  any  given  year  and  using  statistics 
for  the  entire  trench  along  the  Peru-Chile  coast,  a least-squares 
analysis  resulted  in  the  following  expression: 

n(i)  = 0.071+e0'631 

39*  In  using  statistics  for  the  entire  trench  area  along  the 
Peru-Chile  coast,  it  was  assumed  that  the  probability  of  tsunami  occur- 
rence was  uniform  along  the  trench.  This  is  a standard  assumption  for 

17 

earthquake  frequency  analysis.  The  tectonic  justification  of  this 

assumption  lies  in  the  fact  that  a single  sialic  block  or  plate  of  the 
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earth's  crust  or  lithosphere  is  dipping  into  the  Peru-Chile  Trench. 

It  can  be  reasonably  expected  that  the  movement  of  this  single  plate 
is  similar  along  its  entire  length. 

1»0.  In  the  Aleutian  Trench  area,  only  large  tsunamis  occurring 
in  relatively  recent  years  (since  1788)  have  been  recorded  due  to  the 
isolation  of  the  area.  Assuming  an  exponential  coefficient  of  -0.71 
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for  this  trench  area  (determined  in  Reference  18  as  a mean  value  for 
areas  of  the  Pacific  with  the  most  data  on  tsunamis)  and  using  only 
the  reliable  data  for  large  tsunamis  (intensity  greater  than  or  equal 
to  3-5.  Reference  16),  the  following  relation  was  determined  by  a least- 
squares  analysis: 


n(i)  = 0.113e 


-0.71i 


Again t the  probability  of  tsunami  occurrence  was  assumed  to  be  uniform 
along  the  trench. 


Use  of  Numerical  Models 


1*1.  To  relate  the  probability  distribution  of  tsunami  intensities 
to  source  characteristics,  it  was  assumed  that  the  ratio  of  the  source 
uplift  heights  producing  two  tsunamis  of  different  intensities  (as  de- 
fined in  the  previous  section)  is  equal  to  the  ratio  of  the  average 
runup  heights  produced  on  the  coasts  near  these  tsunami  sources.  This 
ratio  is  equal  to  2^1-:'’2^  for  two  tsunamis  with  intensities  i^  and 


1*2.  If  Ha  is  the  wave  height  emitted  in  the  direction  parallel 

to  the  major  axis  of  length  a by  a tsunami  source  with  an  elliptical 

shape  (large  tsunamis  have  historically  had  elliptically  shaped  uplifts) 

and  is  the  wave  height  emitted  in  the  direction  parallel  to  the 

minor  generation  axis  of  length  b , then  experimental  research  of 

tsunami  generation  has  shown  that  ~ a/b  (Reference  22).  For  a 

large  tsunami , H,  can  be  larger  than  H by  a factor  of  as  much  as 
b a 

5 or  6.  Thus,  the  orientation  of  the  tsunami  source  relative  to  the 
area  where  runup  is  to  be  determined  is  very  important;  that  is,  the 
runup  at  a distant  site  due  to  the  generation  of  a tsunami  at  one  loca- 
tion along  a trench  cannot  be  considered  as  being  representative  of  all 
possible  placements  of  the  tsunami  source  in  the  entire  trench  region. 
Hence,  the  Aleutian  and  Peru-Chile  Trenches  had  to  be  segmented  and 
runup  along  the  west  coast  of  the  United  States  determined  for  tsunami 


sources  located  at  the  center  of  each  of  the  segments. 

**3.  The  spatial  size  of  a tsunami  source  was  standardized  because 

there  is  not  an  apparent  correlation  between  tsunami  intensity  and 

spatial  size  of  a tsunami  source.  For  example,  the  I9U6  Aleutian 

tsunami  had  an  uplift  region  of  very  small  spatial  extent,  whereas  the 

1957  Aleutian  tsunami  had  an  uplift  region  that  covered  perhaps  the 
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greatest  spatial  extent  of  any  known  earthquake;  yet  the  19^6  tsunami 
had  the  greater  intensity,  producing,  in  general,  greater  runup  eleva- 
tions in  the  near  and  distant  regions. 

1*1*.  A description  of  the  standard  source  employed,  given  in  Ref- 
erence 1,  represents  a large  tsunami  with  intensity  1*  on  the  modified 
Imamura-Iida  scale.  Certainly,  tsunamis  of  low  intensity  may  have 
smaller  spatial  extents;  however,  large  tsunamis  pose  the  greatest 
threat  to  a distant  area  such  as  the  west  coast  of  the  United  States. 
These  large  tsunamis  can  be  expected  to  have  similar  spatial  extents, 
with  any  spatial  differences  being  unimportant  in  the  far  field  com- 
pared with  the  effects  of  source  orientation  and  vertical  uplift. 

1*5.  Figures  7 and  8 show  the  Aleutian  Trench  divided  into  12 
segments  and  the  Peru-Chile  Trench  into  3 segments.  The  segments  in 
the  Aleutian  Trench  were  approximately  one-quarter  the  length  of  the 
major  axis  of  the  standard  source,  whereas  the  segments  in  the  Peru- 
Chile  Trench  were  approximately  the  length  of  the  major  axis  of  the 
standard  source.  The  standard  source  was  centered  in  each  segment  such 
that  the  major  axis  of  the  source  was  parallel  to  the  trench  axis.  Up- 
lift regions  historically  have  had  such  an  orientation  relative  to 
trench  systems.  The  Aleutian  Trench  was  segmented  much  finer  than  the 
Peru-Chile  Trench  because  the  Aleutian  Trench  is  oriented  relative  to 

the  west  coast  such  that  elevations  produced  on  the  west  coast  are  very 
21* 

sensitive  to  the  location  of  a source  along  the  Trench.  Uplifts 
along  the  Peru-Chile  Trench  do  not  radiate  energy  directly  toward  the 
west  coast  regardless  of  their  position  along  the  Trench.  The  Peru 
and  Chile  sections  of  the  Peru-Chile  Trench  have  constant  orientations 
relative  to  the  west  coast  of  the  United  States;  therefore  elevations 
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Figure  7.  Idealized  axis  of  Aleutian  Trench  showing  12  segments 
and  perimeter  of  standard  uplift  area 
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on  the  west  coast  are  not  very  sensitive  to  source  location  within 
these  sections. 

46.  In  each  of  the  segments  of  the  Aleutian  and  Peru-Chile 
Trenches,  tsunamis  with  intensities  from  2 to  5 in  increments  of  0.5 
were  generated  and  propagated  across  the  deep  ocean  using  the  deep- 
ocean  numerical  model  discussed  in  an  earlier  section.  Tsunamis  with 
intensity  less  than  2 are  too  small  to  produce  significant  runup  on  the 
west  coast.  An  upper  limit  of  5 was  chosen  because  the  largest  tsunami 
intensity  ever  reported  was  less  than  5.'^  Reference  17  indicates  that 
there  is  an  upper  limit  to  the  strain  which  can  be  supported  by  rock 
before  fracture.  Thus,  earthquakes  only  reach  certain  maximum  magni- 
tudes and  tsunamis  can  be  expected  to  have  similar  upper  limits  of  in- 
26  27 

tensity.  Perkins  and  McGarr  have  demonstrated  that  future  earth- 
quakes cannot  have  seismic  moments  (measure  of  earthquake  magnitude  for 
large  earthquakes)  much  larger  than  those  of  earthquakes  that  occurred 
in  recorded  history. 
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Figure  8.  Idealized  axis  of  Peru-Chile 
Trench  showing  3 segments 


U7.  The  wave  forms  propagated  to  the  west  coast  by  the  deep- 
ocean  propagation  numerical  model  were  used  as  input  to  the  nearshore 
propagation  numerical  model.  Each  wave  form  was  propagated  from  a 
water  depth  of  500  m to  shore  using  the  nearshore  model  and  one  of  the 
grids  shown  in  Figure  3.  Thus,  at  each  grid  location  on  the  shoreline 
of  the  west  coast,  there  was  a group  of  105  wave  forms — seven  wave 
forms  (for  intensities  from  2 to  5 in  increments  of  0.5)  for  each 


segment  of  the  Aleutian  and  Peru-Chile  Trenches.  Each  of  t.he*e  wave 
forms  had  an  associated  probability  equal  to  the  probability  that  a 
certain  intensity  tsunami  will  be  generated  ir.  a particular  segment  of 
a trench  region. 


Effect  of  Astronomical  Tides 


U8.  The  maximum  "still  water"  elevation  produced  during  tsunami 
activity  is  the  result  of  a superposition  of  tsunami  and  tidal  wave 
forms.  Therefore,  th°  statistical  effect  of  astronomical  tides  on 
total  tsunami  runup  must  be  included  in  the  predictive  scheme  presented 
in  this  report.  Since  the  wave  forms  calculated  by  the  nearshore  model 
did  not  have  a simple  form  (e.g. , sinusoidal),  the  statistical  effect 
of  the  astronomical  tide  on  tsunami  runup  had  to  be  determined  through 
a numerical  approach. 

1*9.  The  wave  forms  calculated  by  the  nearshore  numerical  model 
extended  ever  a period  of  time  of  approximately  2 hr.  Three  or  four 
wave  crests  (the  largest  waves  in  the  tsunami)  arrived  during  this  time. 
Smaller  waves  arriving  at  later  times,  however,  have  often  persisted  for 
days  during  historical  tsunamis.  An  analysis  of  tide  gage  records  of 
the  i960  and  l°6k  tsunamis  on  the  west  coast  indicated  that  these 
smaller  waves  have  amplitudes  on  the  average  of  Uo  percent  of  the  maxi- 
mum wave  amplitude  of  the  tsunami;  therefore  a sinusoidal  group  of  these 
small  waves  was  e.dded  to  each  of  the  calculated  wa.ve  forms  so  that  the 
total  wave  form  extended  over  a 2U-hr  period.  These  sma.ller  waves  are 
important  for  locations  where  tsunami  waves  are  fairly  small  compared 
with  tidal  variations.  At  such  locations,  the  maximum  combined  tsunami 
and  astronomical  elevation  occurs  during  the  maximum  tidal  elevation. 

50.  A computer  program  was  developed  to  predict  time-histories 

of  the  astronomical  tides  at  all  grid  locations  on  the  west  coast  of 

the  United  States.  The  program  was  based  upon  the  harmonic  analysis 

methods  used  in  the  past  by  the  National  Ocean  Survey  for  mechanical 
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tide-predicting  machines.  Tidal  constants  available  from  the  National 
Ocean  Survey  were  used  as  input  to  the  computer  program.  A year  of 


tidal  elevations  was  then  predicted  for  grid  locations  all  along  the 
west  coast.  The  year  19 6k  was  selected  because  all  the  major  tidal 
components  had  a node  factor  of  approximately  1.00  during  this  year, 
thus  making  it  an  average  year.  The  node  factor  is  associated  with  the 
revolution  of  the  moon's  node  and  has  an  18. 6-year  cycle.  Since  a 
tsunami  can  arrive  at  any  time  during  this  18. 6-year  period  (arrival  at 
a low  of  the  node  factor  is  equally  likely  as  an  arrival  at  a high), 
the  statistical  effect  of  the  varying  node  factor  is  small  and  an  aver- 
age value  should  be  used.  The  statistical  effect  of  the  varying  node 
factor  on  the  predicted  runup  elevations  can  be  shown  to  be  a small 
fraction  of  an  inch  using  the  approach  discussed  in  paragraphs  56  and 
57  (with  a a for  the  nodel  variation  approximately  equal  to  0.5  in. 
squared) . 

51.  The  tidal  time-histories  calculated  at  each  of  the  nearshore 
numerical  model  grid  points  along  the  west  coast  were  then  subdivided 
into  15-min  segments.  Each  of  the  105  2U-hr  wave  forms  was  allowed  to 
arrive  at  the  beginning  of  each  of  these  15-min  segments  and  then  super- 
imposed upon  the  astronomical  tide  for  the  2k-hr  period.  The  maximum 
combined  tsunami  and  astronomical  tide  elevation  over  the  2U-hr  period 
was  determined  for  tsunami  wave  forms  arriving  during  each  of  these  15- 
min  starting  times.  Each  of  these  maximum  elevations  had  an  associated 
probability  equal  to  the  probability  that  a certain  intensity  tsunami 
would  be  generated  in  a particular  segment  of  the  two  trench  regions 
and  arrive  during  a particular  15-min  period  of  a year. 

52.  The  many  maximum  elevations  with  associated  probabilities 
were  used  to  determine  cumulative  probability  distributions  of  combined 
tsunami  and  astronomical  tide  elevations.  The  maximum  elevations  were 
ordered  and  probabilities  summed,  starting  with  the  largest  elevations, 
until  a desired  probability  was  obtained.  The  elevation  encountered 
when  the  summed  probabilities  reach  a desired  value  P was  the  eleva- 
tion that  is  equaled  or  exceeded  with  an  average  frequency  of  once 
every  1/P  years.  Thus,  when  the  summed  probabilities  reached  the  value 
0.01,  the  elevation  associated  with  the  last  probability  summed  was  the 
100-year  elevation. 
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53.  100-  and  500-year  elevations  were  determined  at  all  grid 

points  of  the  nearshore  numerical  grids.  Smooth  curves  were  then  used 
to  connect  all  these  discrete  elevations,  resulting  in  the  continuous 
elevation  curves  shown  in  Plates  1-30. 

Comparison  with  Local  Observation  Predictions 

5**.  Crescent  City,  California,  was  the  only  location  on  the  por- 
tion of  the  west  coast  of  the  United  States  considered  in  this  report 
(that  is,  the  west  coast  minus  Southern  California,  Monterey  Bay,  San 
Francisco  Bay,  and  Puget  Sound)  which  had  sufficient  historical  data  of 
tsunami  activity  to  allow  frequency  of  occurrence  predictions  based 
upon  local  historical  observations  (note,  however,  discussion  for  San 
Francisco  in  next  section).  Wiegel'*'^  made  such  predictions  of  tsunami 
height  (trough  to  crest  height)  for  the  period  from  1900  to  1965-  He 
predicted  100-  and  500-year  heights  at  Crescent  City  of  approximately 
25-6  ft  and  U3.2  ft  (extending  his  curve)  respectively.  If  the  crest 
amplitude  is  taken  to  be  one-half  the  total  height,  the  100-  and  500- 
year  elevations  are  12.8  ft  and  21.6  ft,  respectively.  However,  the 
crest  amplitude  at  Crescent  City  is  typically  greater  than  one-half  the 
wave  height  (e.g.,  the  crest  amplitude  of  the  largest  wave  of  the  196U 
tsunami  was  approximately  60  percent  of  the  total  height).  If  Wiegel's 
analysis  is  applied  to  historical  crest  elevations  instead  of  heights, 
the  100-  and  500-year  elevations  are  found  to  be  15.1*  ft  and  26. U ft, 
respectively.  Furthermore,  the  analysis  now  can  be  applied  to  the 
longer  time  period  from  1900  through  1977.  The  100-  and  500-year  ele- 
vations based  upon  this  longer  time  span  are  lU.5  ft  and  25.5  ft, 
respectively. 

55«  Plate  17  shows  100-  and  500-year  elevations  at  Crescent  City 
of  13.1  ft  and  2U.9  ft,  respectively.  These  values  compare  very  favor- 
ably with  the  lU.5-  and  25. 5- ft  elevations  determined  from  historical 
data  for  the  period  of  time  from  1900  through  1977.  The  elevations  pre- 
dicted by  the  analysis  based  upon  the  local  historical  data  are  probably 
somewhat  larger  than  the  elevations  predicted  in  this  report  because  the 
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short  time  period  (relative  to  100  and  500  years)  from  1900  through  1977 
includes  the  exceptionally  active  years  of  tsunami  generation  from  19^6 
through  1961k  The  three  greatest  intensity  tsunamis  generated  in  the 
Aleutian-Alaskan  region  since  at  least  1788  occurred  during  this  time 
span.  The  i960  Chilean  tsunami  was  the  greatest  intensity  tsunami 
generated  in  the  Peru-Chile  region  since  at  least  1562. 

56.  The  13. 1-  and  2l4.9-ft  elevations  predicted  for  Crescent  City 
using  the  techniques  discussed  in  this  report  are  not  quite  comparable 
to  the  elevations  predicted  using  Wiegel's  analysis  based  upon  local 
historical  data  because  the  effect  of  the  astronomical  tide  has  been 
included  in  the  elevations  predicted  in  this  report  but  not  in  Wiegel's 
analysis.  However,  the  statistical  effect  of  the  astronomical  tide  on 
the  tidal  elevation  is  small  at  Crescent  City  because  tsunamis  are  so 
large.  Reference  1 shows  that  statistical  effect  of  the  astronomical 
tide  for  a location  where  tsunami  waves  are  large  (so  that  one  wave  is 
much  larger  than  any  other)  is  to  increase  the  predicted  elevation  by 
an  amount  equal  to 
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where  a is  the  tidal  variance  and  equals 


CO 


l 

m=l 


m 


where  C is  equal  to  the  tidal  constituent, 

m 

57.  The  term  a is  given  by  the  following  expression: 

P(g)  = Ae-a& 

where 

A = constant  coefficient 

2 = the  runup  at  any  time  above  local  mean  sea  level 

P(S)  = the  cumulative  probability  distribution  for  runup  at  a 
given  site  being  equal  to  or  exceeding  S due  only  to 
the  maximum  wave  of  the  tsunami. 


This  form  of  the  probability  distribution  is  well  established  for 
tsunamis  (see  paragraphs  35-Uo). 

58.  Based  upon  the  tidal  constituents  predicted  by  the  National 
Ocean  Survey,  o = 7.1  for  Crescent  City.  Using  Wiegel's  data,  the 
period  of  time  from  1900  through  1977  gives  a = O.1U5  . Therefore, 
the  astronomical  tide  contributes  approximately  0.5  ft  to  the  elevations 
predicted  for  Crescent  City. 

San  Francisco  Elevation  Predictions 

59-  Frequency  of  occurrence  calculations  for  San  Francisco  Bay 

2 

were  made  in  an  earlier  report  ; however,  a large  section  of  San  Fran- 
cisco Bay  was  included  in  the  nearshore  numerical  model  grid  (Figure  3) 
used  in  this  study  in  order  that  the  effect  of  the  bay  on  elevations 
outside  the  bay  could  be  properly  simulated.  Thus,  elevation  calcula- 
tions were  made  in  this  study  at  the  location  of  the  Presidio  tide  gage 
in  San  Francisco  Bay.  Plate  9 shows  the  100-  and  500-year  elevation 
predictions  at  this  location  to  be  6.1  ft  and  10.0  ft,  respectively. 
Reference  2 predicted  elevations  of  8.0  and  15.5  ft  near  the  Presidio 
tide  gage.  However,  the  simple  one-dimensional  nearshore  solution 
employed  in  Reference  2 could  only  propagate  the  leading  wave  of  the 
tsunami  to  a location  near  the  mouth  of  San  Francisco  Bay  and  then  it 
had  to  be  assumed  that  the  elevation  at  the  Presidio  tide  gage  was  the 
same  as  this  calculated  elevation.  Thus,  the  elevations  predicted  for 
the  tide  gage  in  Reference  2 really  are  elevations  for  San  Francisco 
just  outside  the  bay;  consequently,  the  elevations  predicted  in  Refer- 
ence 2 at  the  Presidio  tide  gage  and  throughout  the  bay  were  too  large. 
However,  the  percentage  error  in  elevations  for  the  interior  of  San 
Francisco  Bay  will  be  less  than  the  percentage  error  for  the  Presidio 
tide  gage,  since  tsunami  elevations  are  smaller  in  the  bay  and,  conse- 
quently, the  effect  of  the  astronomical  tides  proportionately  greater. 

In  fact,  for  much  of  San  Francisco  Bay  the  error  in  the  tsunami  eleva- 
tion would  probably  be  so  small  compared  with  the  contribution  from  the 
astronomical  tide  that  it  would  be  negligible.  If  the  reduction  of  a 
predicted  elevation  in  San  Francisco  Bay  by  a percentage  equal  to  the 


percentage  error  at  the  Presidio  tide  gage  would  have  a significant  im- 
pact on  a development,  then  recomputations  based  upon  the  results  of 
this  study  might  be  warranted. 

60.  There  are  sufficient  historical  data  at  the  Presidio  tide 
gage  in  San  Francisco  to  make  elevation  predictions  based  upon  local 
historical  data.  Using  a logarithmic  distribution  and  Wiegel's  wave- 
height  data  for  San  Francisco  over  the  period  1900  through  1977  yields 
100-  and  500-  year  heights  (trough  to  crest)  of  8.0  ft  and  lb. 2 ft, 
respectively.  If  the  crest  amplitude  is  one-half  the  height,  this 
yields  100-  and  500-year  elevations  of  lt.0  ft  and  7.6  ft,  respectively. 

A similar  analysis  just  based  upon  historical  crest  amplitudes  (with 
no  assumption  concerning  the  crest  to  trough  ratio)  yields  100-  and 
500-year  elevations  of  lt.lt  ft  and  7.8  ft,  respectively. 

61.  It  is  difficult  to  compare  the  elevations  predicted  using 
local  historical  data  with  the  elevations  of  6.1  and  10.0  ft  calculated 
in  this  report,  because  the  effect  of  the  astronomical  tides  (contained 
in  the  calculations  of  this  report)  cannot  easily  be  estimated,  as  was 
the  case  for  Crescent  City.  Tsunamis  recorded  at  the  Presidio  tide  gage 
in  San  Francisco  are  known  to  persist  at  fairly  substantial  levels  for 
extended  periods  of  time,  apparently  due  to  some  oscillation  phenomenon. 
For  example,  during  the  1961t  tsunami  there  were  six  waves  of  approxi- 
mately equal  or  greater  amplitude  than  the  initial  wave  (in  addition 

to  smaller  waves)  during  the  first  8 hr  of  the  tsunami.  Oscillations 
approximately  bo  percent  of  the  amplitude  of  the  initial  wave  persisted 
for  at  least  another  2b  hr.  The  tide  range  at  San  Francisco  is  5.7  ft 
with  mean  higher  high  water  (mhhw)  2.7  ft  above  msl  and  mean  high 
water  (mhw)  1.0  ft  above  msl.  At  least  one  of  the  main  waves  of  a 
tsunami  at  San  Francisco  probably  superimposes  upon  either  mhw  (with 
the  smaller  oscillations  then  adding  on  to  mhhw)  or  mhhw.  Thus,  the 
effect  of  the  tides  is  probably  to  contribute  1 to  2.7  ft  to  the  total 
elevation.  The  average  of  these  tidal  elevations  is  1.85  ft  and  if 
this  elevation  is  added  to  the  I4.lt-  and  7-8-ft  elevations  based  upon 
local  historical  data,  the  resulting  elevations  of  6.25  and  9.65  ft  are 
similar  to  the  6.1-  and  10.0-ft  elevations  calculated  in  this  report. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

62.  The  deep-ocean  and  nearshore  numerical  models  used  in  this 
report  accurately  simulated  tsunami  generation,  deep-ocean  and  nearshore 
propagation,  and  interaction  with  coastlines.  The  nearshore  numerical 
model  was  superior  to  the  analytical  techniques  used  for  nearshore 
propagation  in  previous  studies  (References  1 and  2)  at  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  because  it  was  two- 
dimensional,  solved  equations  which  included  nonlinear  advective  terms 
and  bottom  stress  terms,  and  handled  arbitrary  incident  wave  forms.  It 
may  be  advisable  to  use  the  nearshore  numerical  model  to  recompute  ele- 
vations for  the  areas  considered  in  References  1 and  2.  However,  since 
the  elevations  computed  in  these  earlier  studies  were  in  general  not 
very  large  (because  the  areas  considered  were  protected  from  tsunamis), 
there  may  not  be  a need  for  the  refinement  in  nearshore  calculations 
provided  by  the  nearshore  numerical  model.  Nonlinear  and  bottom  stress 
effects  also  are  less  important  for  smaller  tsunamis.  Thus,  neglect  of 
these  effects  in  References  1 and  2 probably  was  not  generally 
significant. 

63.  The  techniques  (which  are  based  upon  historical  data  from 
source  regions  and  numerical  model  calculations)  used  in  this  report 
to  predict  tsunami  runup  provided  good  estimates  of  100-  and  500-year 
runup  levels  for  the  west  coast  of  the  continental  United  States . The 
100-  and  500-year  predictions  at  Crescent  City  and  San  Francisco,  Cali- 
fornia, based  upon  the  techniques  presented  in  this  report,  were  shown 
to  agree  quite  well  with  predictions  based  upon  historical  observations 
of  tsunami  activity  at  these  two  locations.  These  were  the  only  loca- 
tions along  the  coastline  considered  in  this  report  which  had  sufficient 
local  historical  data  to  allow  predictions  based  upon  local  observations. 

Recommendations 

61*.  The  runup  elevations  predicted  in  this  report  can  be  assumed 
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to  equal  the  tsunami  shoreline  elevation  for  most  of  the  shoreline  pre- 
sented in  Plates  1-30.  The  similarity  between  runup  elevations  and 
elevations  at  the  shoreline  is  illustrated  in  the  detailed  surveys  of 

tsunami  inundation  that  have  been  performed  following  historical 

14 

tsunamis.  For  example,  Magoon  reported  flooding  to  about  the  20- ft 

contour  above  mllw  and  elevations  at  the  shoreline  of  about  20  ft  (mllw) 

for  the  1964  tsunami  at  Crescent  City,  California.  Wilson  and  Torum^ 

reported  that  the  20- ft  (mllw)  runup  at  Valdez,  Alaska,  for  the  1964 

tsunami  checked  "well  for  consistency  with  water  level  measurements 

made  on  numerous  buildings  throughout  the  town."  Similar  comments  were 

made  by  Brown‘d  in  reference  to  survey  measurements  of  30-ft  (mllw) 

runup  at  Seward,  Alaska,  for  the  1964  tsunami.  Runup  and  elevation  at 

the  shoreline  were  similar  at  nine  locations  in  Japan  as  recorded  by 
30 

Nasu  in  surveys  following  the  1933  Sanriku  tsumani.  Runup  and  eleva- 
tion at  the  shoreline  were  also  similar  for  the  borelike  waves  at  Hilo, 

31 

Hawaii,  during  the  i960  tsumani. 

65.  There  are  locations  where  time-dependent  effects  (for  example, 

lack  of  sufficient  time  to  completely  flood  extensive  low-lying  or 

estuarine  areas)  or  two-dimensional  effects  (for  example,  flow 

divergence  or  convergence)  cause  tsunami  runup  elevations  not  to  be 

equal  to  elevations  at  the  shoreline.  It  is  recommended  that  for  these 

areas  inundation  limits  be  determined  using  a numerical  model  developed 
32 

recently  at  WES.  This  model  is  capable  of  handling  land  flooding  for 
bays,  harbors,  developed  areas  such  as  cities,  large  low-lying  areas, 
sand-dune  protected  areas,  and  other  areas  where  there  are  topographi- 
cal, roughness,  or  coastline  variations. 
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APPENDIX  A:  NOTATION 

a Length  of  major  axis  of  elliptical  source,  miles 
A Probability  constant 

b Length  of  minor  axis  of  elliptical  source,  miles 

1/2 

C Chezy  coefficient,  ft  /sec 


of  ellipse,  ft 
of  ellipse,  ft 

H Average  runup  over  a coast,  m 

avg 

1 Tsunami  intensity 

n(  ) Tsunami  probability  function 

P(  ) Cumulative  probability  distribution  for  runup 
t Time,  sec 

x 2 2 

T Bottom  stress  in  x-direction,  ft  /sec 
B 2 2 

T y Bottom  stress  y-direction,  ft  /sec 
B 

U Water  velocity  in  x-direction,  ft/sec 
V Water  velocity  in  y-direction,  ft /sec 
x Spatial  coordinate,  ft 

y Spatial  coordinate,  ft 

3 Water-surface  elevation  above  mean  sea  level  (msl),  ft 
a Exponential  constant,  ft  ^ 
n Wave  height,  ft 

o Model  variation 

2 2 
a Tidal  variance,  ft 


m 

e 

g 

h 


th 


Amplitude  of  m tidal  constituent,  ft 

2 

Acceleration  due  to  gravity,  ft/sec 
Water  depth , ft 

Wave  height  in  direction  of  major  axis 
Wave  height  in  direction  of  minor  axis 


A1 


